CERN-TH/97-370 
On the Index Theorem for Wilson Fermions 



P. Hernandez 

Theory Division, CERN, 1211 Geneva 23, Switzerland. 

Abstract 

We consider a Wilson-Dirac lattice operator with improved chiral properties. We 
show that, for arbitrarily rough gauge fields, it satisfies the index theorem if we identify 
the zero modes with the small real eigenvalues of the fermion operator and use the 
geometrical definition of topological charge. This is also confirmed in a numerical study 
of the quenched Schwinger model. These results suggest that integer definitions of the 
topological charge based on counting real modes of the Wilson operator are equivalent 
to the geometrical definition. The problem of exceptional configurations and the sign 
problem in simulations with an odd number of dynamical Wilson fermions are briefly 
discussed. 
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1 Introduction 



The connection between gauge field topology and fermion zero modes is expected to 
have important physical consequences on the non-perturbative dynamics of baryon 
number violation in the SM or the breaking of the singlet chiral symmetry in QCD. 
The study of these effects, however, requires non-perturbative techniques and one would 
expect that Monte Carlo methods on the lattice would ultimately be best suited to 
it. Unfortunately there is no proof of the Atiyah-Singer index theorem on the lattice, 
in spite of which, a big effort has been devoted to the measurement of the topological 
susceptibility ffl]. The main motivation for this measurement in SU(3) is the relation of 
this quantity to the rj mass, through the continuum formula of Witten and Veneziano 
0| . This formula relies on the validity of the index theorem, so it is important to 
understand to what extent this theorem survives on the lattice, where the connection 
between topology and fermion zero modes is not clear. 

It has been known for some time that there are remnants of the index theorem on the 
lattice for Wilson fermions. In || it was observed that in smooth gauge configurations, 
in which the geometrical definition of topological charge pi] takes an integer value 
Q 9eo , the Wilson-Dirac operator has Q geo small and exactly real eigenvalues with the 
appropiate chirality. The small real eigenvalues seem to play the role of the continuum 
zero modes. We will refer to the connection between the small real eigenvalues of 
the Wilson-Dirac operator and the geometrical definition of topological charge as the 
"lattice index theorem" (LITh) . Unfortunately, when rough gauge fields are considered 
this connection seems to be lost. 

It is well known that the geometrical definition of topological charge in four- 
dimensional Yang- Mills theories has dislocations. They are small , 0(a), objects, which 
carry topological charge, but have such a small plaquette action that they destroy the 
proper scaling of topological quantities constructed out of the geometrical charge. It 
has been argued that the geometrical definition of topological charge, being constructed 
as an integral of a local density, is more sensitive to fluctuations of the order of the 
lattice spacing than the number of fermion "zero" modes, which are non-local ||. For 
this reason, there have been several proposals Q H to measure topology by measuring 
the chiral charge of the small eigenvalues of the Wilson-Dirac operator. The hope is 
that these definitions are free of dislocations and topological quantities show the ex- 
pected scaling. If this were the case, then it is necessary that the LITh is violated, 
even in an average sense, as we approach the continuum limit, because the fermionic 
charge should scale differently than the geometrical one. The experimental fact is that 
the LITh is violated for rough gauge fields, but satisfied for smooth ones. This does 
not allow us to draw any conclusion a priori on whether the LITh will be satisfied in 
the continuum limit (in an average sense), because rough fields are important even in 
this limit, since they are responsible for renormalization. 

In this letter, we present a study of the LITh for a new fermionic action that was 
originally proposed to deal with chiral gauge theories ]?[], which has improved chiral 
properties at any fixed lattice spacing. We apply it here to the study of the Schwinger 
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model. The new action is constructed by interpolating the gauge variables smoothly, 
gauge invariantly, and locally to a finer lattice ||, in which the fermion determinant 
and propagator are defined in terms of the standard Wilson-Dirac operator. It can 
be shown that the factor of fine-graining controls the violations of chirality Q like a 
power, to all orders in the gauge coupling (in the appendix, we present an explicit one- 
loop calculation in a £7(1) model in four dimensions of the additive renormalization to 
the fermion mass, which shows the expected suppression.). The improvement in the 
validity of LITh is however much more dramatic than this power suppression might 
indicate. We find that, with a simple factor of fine-graining of 1/2, the violations of 
LITh are absent for arbitrarily large couplings. This strongly indicates that measuring 
topological charge by looking at the real eigenvalues of the Wilson operator (when it 
has been properly improved) is equivalent to the geometrical definition. 

If this is confirmed in four dimensions, the problem of dislocations should be handled 
by improving the plaquette action |ll], [l2|, |l3| rather than by using a different integer 
definition of the topological charge. Any such definition is not sensitive to topological 
objects of size roughly of 0(a), since a is the only cutoff scale in the problem. Different 
definitions might have a slightly different cutoff, but generically there is no reason to 
expect that some of them will be affected by dislocations and not others, as long as 
their cutoff is of the same order []. The improved gauge action described in [13] is the 
natural choice to combine with the improved fermionic action described here in order 
to ensure the proper scaling of the topological susceptibility constructed out of the 
geometrical or fermionic charges. 

In section 2, we review some known facts about the connection between zero modes 
in the continuum and real modes of the Wilson-Dirac operator and show that a lattice 
index theorem can be defined, which holds for every gauge configuration if we use the 
improved Wilson-Dirac operator. In section 3, we present results on the LITh for the 
standard action. In section 4, the improved action is described and the results on 
the LITh are presented. The small real eigenvalues of the Wilson-Dirac operator are 
related to the so-called exceptional configurations and are also responsible for the sign 
flips in the determinant of this operator, which make dynamical simulations with an 
odd number of fermions problematic. In section 5, we discuss the relevance of our 
results to these problems and conclude. 



2 Lattice Index Theorem 

We first review the known arguments that suggest the identification of the zero modes 
of the continuum Dirac operator with the exactly real eigenvalues of the Wilson-Dirac 
operator on the lattice §, g, |, [lj, p]. 

Continuum zero modes have a well-defined chirality. As was shown in f^, f|, the 
eigenmodes of the lattice Wilson-Dirac operator with a non-vanishing chirality are 



^The fermionic charge defined in jj| is not an integer so this argument does not apply to it. 
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necessarily real. This is easy to see by realizing that the Wilson-Dirac operator satisfies, 

75 ^75 =P\ (1) 

which implies that eigenvalues come in complex conjugate pairs, i.e. if V{ is an eigen- 
vector on the right of p with eigenvalue Aj, then vjjs is an eigenvector on the left with 
eigenvalue A*. Then, it follows that, if vj^Vi ^ 0, 

A* = A* (2) 

must be satisfied and, vice versa, if Aj is complex then vjj^Vi = 0. Generically, for real 
eigenvalues vj^Vi = 0(1). In fact as found in ||, this value is very close to ±1 for 
smooth backgrounds. 

Another indication that real eigenvalues might have a topological origin is their 
stability under perturbations. Consider the Wilson-Dirac operator P =f)(°) + e P^\ 
Both terms satisfy the property ([l]). Let Vi be the eigenvectors of P^ and v[ the 
perturbed ones. The chirality of the perturbed eigenvectors to leading order in e is 

v'hbv'i = vjj 5 Vi + 0(e) 2 . (3) 

If Vi is an eigenvector ofi^ ) with real eigenvalue, the first term on the right-hand side 
is of O(l). Then a very large perturbation would be needed for v[ to correspond to a 
complex eigenvalue oip since, in this case, the left-hand side of (|3|) vanishes. Similarly, 
for a complex eigenvalue of P^ with eigenvector Vi to become real, the left-hand 
side of (||) should be of 0(1), while the first term on the right-hand side vanishes. A 
perturbation of 0(1) is needed. 

The continuum index theorem implies that the number of zero modes with positive 
chirality minus the number of eigenvalues with negative chirality should be equal to the 
topological charge. As the lattice definition of the topological charge we will choose the 
geometrical definition of Liischer Q, Q geo . This charge is defined as the naive charge 
of a continuum gauge configuration obtained by smoothly interpolating the lattice 
configuration to the continuum. A continuum gauge field a^(x) can be constructed out 
of the link variables, which satisfies 

U^s) = e l $s + * dxa » {x \ (4) 

and which transforms covariantly under a lattice gauge transformation. For Yang-Mills 
in four dimensions, Q geo is then defined as 

Qgeo — _____ J aftxfpvfpv, (5) 

where f^ u is the field strength of the field a M . 

Numerically, it has been found in all previous investigations (see for example ||) 
that the number of real eigenvalues of the Wilson-Dirac operator is a multiple of 2 d and 
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their net chirality vanishes. As we will see, this is the consequence of fermion doubling 
in non-trivial backgrounds. On the other hand, the Wilson term is responsible for 
giving a large mass to the doublers, which translates in the fact that the 2 d sets of real 
eigenvalues cluster in d + 1 regions of the real axis. For smooth backgrounds, the real 
eigenvalues appear near the points A = 2rnj/a, where r is the Wilson coupling, a is 
the lattice spacing, no = for the physical modes, and n,, for i = 1, 2 d — 1, counts 
the number of momentum components that are equal to ir in the ith doubler corner of 
the Brillouin zone. We then define the physical region as the interval of the real axis 
S p = [0, r/a], where we expect to find the modes corresponding to the the physical 
fermion, and the doubler regions as S^. = [(2nj — l)r/a, (2n« + \)r/d\. Notice that the 
regions corresponding to different doublers with the same nj coincide. 

This pattern of the distribution of the real eigenvalues of D is actually easy to prove 
for smooth backgrounds (to my knowledge this has not been proved before). Let us 
consider the lattice Wilson-Dirac operator (without bare mass): 

p = lE^^+^-iE^;* ( 6 ) 

H fj, 

where the covariant are given by D+ = 5 s r s+ pUfa) - 5 S , S , D~ = 6 S > S - Ufa - /t)<Va-j«- 
We also define the continuum Euclidean operator: 

where is the continuum field that satisfies eq. (||) and whose topological charge 
is Q geo . The continuum index theorem ensures that TJ) C has zero modes with a net 
chirality Q geo . On the other hand, the lattice operator (^) can be expanded (at small 
lattice momentum) in the lattice spacing, a, 

0=0 C + a + (8) 

The terms of 0(a) are small, because we are considering smooth backgrounds and small 
lattice momentum. Then the zero modes of p c become small real eigenvalues of p. 
This is because the perturbation of the small 0(a) corrections can only move the zero 
modes along the real axis, but not make them complex, as we have explained before. 

Considering the doubler fermions, we must perform a similar expansion in a, but 
around the appropiate corner of the Brillouin zone. This can be achieved by first 
performing a unitary transformation of the lattice operator and then expanding naively 
in a, 

= Vl $>V di = r + a DID* + (9) 

where 

v d,( s ) = II(^r o + ^75^71-) exp(i Kis), (10) 

V 
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with iii = 1, ..,d. Ki is the ith doubler momentum (e.g. the lightest doublers in 
four dimensions correspond to Ki = (ir, 0, 0, 0), (0, ir, 0, 0), (0, 0, ir, 0), (0, 0, 0, ir) with 
Hi = 1). It is easy to check that applying the Vi rotation to the free lattice operator, 
we map the small momentum region into the region surrounding Kj. Again, for smooth 
backgrounds, the 0(a) terms are small, so the zero modes of p c become real modes 
of with eigenvalues 2rrii/a + O(a) (again the shift of the eigenvalues by the O(a) 
effects is along the real axis). On the other hand, applying to the eigenvectors of 
we obtain eigenvectors of the 1J) with the same eigenvalues, which implies that jp 
also has real eigenvalues at 2rrii/a + 0(a). This shows that for each zero mode of JD C 
in a smooth background, there are 2 d real modes of Tf). Furthermore, if we define the 
lattice chirality of each eigenvector Vi as 

Xi = sign(vjj 5 Vi). (11) 

the net chirality in each of the 2 d sets is (— l) ni Q geo , since the Vi for odd rii are 
chirality- flipping matrices. 

The lattice version of the index theorem then becomes, 

N R - N L = Q9 e °, (12) 

where N^l are the real eigenvalues in S p with positive and negative chirality. 

As we have shown, ( |l2|) is satisfied for smooth backgrounds, but not necessarily 
when rough gauge fields are considered. In general rough gauge configurations are 
important even at small coupling, since they are responsible for renormalization. It 
is then not clear whether the two integer definitions of topological charge, Q geo and 
Nr — Nl, and other quantities constructed out of them, are equivalent up to 0(a) 
effects or if they differ even in the continuum limit. 

In two dimensions, ( |l2|) will be exact at small enough coupling. The reason is that 
the effects of the higher dimensional operators of (||) and @ are truly O(a), because 
the integration over gauge fields in two dimensions can at most give logarithmically 
divergent contributions, which cannot compensate the naive power suppresion. This 
is of course due to the superrenormalizability of the theory. In four dimensions, this 
is not so obvious. The higher dimensional operators in (||) and @, upon gauge field 
integration, will induce for example a divergent renormalization of the fermion mass, 
which will shift all the zero modes of IJ) C by a constant (in general it can be a differ- 
ent constant with alternating signs in S p and in S^)- Let's call the averages of the 
eigenvalues in S p and S^, m c and respectively. For smooth backgrounds, we saw 
that m c = and = 2rrii/a. Doubler decoupling requires that this hierarchy be 
maintained under renormalization, in such a way that the physical mass of the lightest 
doubler Ai = — ra c ~ 1/a. On the other hand, in order to ensure the separation 
between the two lightest sectors, the dispersion of the real eigenvalues around their 
means m c and m^, that we call a c and a^, should be much smaller than this. It 



has been conjectured P][1C] that this dispersion goes to zero in the continuum limit. 



Although this seems reasonable, we have not been able to find a rigorous proof for 
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Figure 1: Quenched distribution of the real eigenvalues values of the Wilson-Dirac operator 
in a two-dimensional U(l) model at (5 — 3.0 (solid) and 1.0 (dashed), at fixed physical 
volume. 

it. On the other hand, if it is true, this would seem to imply that the effect of the 
higher dimensional operators in (|8|) and @ is, up to corrections vanishing in the con- 
tinuum limit, a shift of the zero modes of p c to the values m c and m^. (If these 
higher dimensional operators would induce non-vanishing effects other than the mass 
renormalization, this would translate generically in a dispersion of the real eigenvalues 
which would not vanish in the continuum limit.) A consequence of this is that the two 
charges Nr — Nl and Q 9eo are the same in the continuum limit (in particular if one 
has dislocations so does the other). 

Actually, it is quite possible that all the violations of the LITh at strong coupling 
are related to the fact that the gap between S p and 5^ closes. This gap is illustrated in 
Fig. 1, where we show the distribution of real eigenvalues of the standard Wilson-Dirac 
operator resulting from a quenched simulation at (3 = 3.0, 1.0 in a two-dimensional U(l) 
model. At the larger (3 value, the distribution clearly signals the expected regions of 
the real axis. However at (3 = 1.0 the mixing between these regions is very large, and 
as a result assigning real modes to the physical or doubler sectors becomes ambiguous. 
This might be considered as an effect of the non-decoupling of the doublers in the 
measurement of topology. 

If the values of <7 P)( ^ are truly 0(a), it is expected that they could be further reduced 
through standard improvement techniques to make them O(a) 2 , in such a way that the 
gap between the different sectors be ensured for larger couplings. A first investigation 
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of this has been presented in and surprisingly a negative result has been found: the 
Clover term (which is the only operator of the appropriate dimensions and symmetries) 
does not seem to change a p . In this letter, however, we propose a new improvement. 
We use "improved" in a somehow loose sense: the action is "improved" in the sense 
that it has improved chiral properties. Since the 0(a) corrections of the standard 
Wilson action are chirality-breaking ones, the new action has smaller 0(a) corrections. 
As we will see in section 4, this action satisfies 

at small momentum, with e being a small parameter (at any gauge coupling) to be 
defined in section 4. This implies that a v and a^, and also |A — 2r/a\ are controlled 
by 0{e) and not by the lattice spacing, ensuring in four dimensions a clean splitting 



between doubler sectors. Actually, (13) implies automatically that for small enough 
e, the (12) is satisfied for arbitrarily large couplings. The proof is identical to the one 
we gave for the standard action on smooth background gauge fields. Futhermore, one 
can also show that the chirality of the real modes is ±1 + O(e), approaching closer the 
continuum behaviour. 

Before discussing our results, it is worth pointing out that the overlap method 
to measure the topological charge || is also equivalent to the left-hand side of eq. 
(|l2|). The charge in Q is related to the number of level crossings of the Hamiltonian 
H(fi) = 75 (JJ) — /i). Clearly at the values of // at which H(fi) has a vanishing eigenvalue, 
the operator p — fi has a zero mode, with the same eigenvector, since this operator 
and H(/j,) have a common kernel. On the other hand, the eigenvectors of P^ and p 
are the same for any fi, so the zero modes of the first correspond to real eigenvalues of 
the second with eigenvalue fi. Then the level crossings of H(fi) occur exactly at the 
values of fi that coincide with the real eigenvalues of p. The charge is then defined as 
the number of positive chirality crossings, minus the number of negative ones, which 
is exactly the left hand side of eq. (|T2|). Also in this case, one has to define what the 
physical region is and look at the crossings that occur only at /Li < M, where M is a 
large enough scale to separate the physical from the doubler sectors ( in Q it is taken 
as M = 1/a, for r = 1, which coincides with S p ). Clearly, the previous reasoning is 
then also applicable to the overlap method. In particular the problem of the absence 
of a gap between the real modes of the physical and doubler sectors will be the same. 

3 LITh with the standard action 

We have considered a quenched two-dimensional compact U(l) model. The gauge 
action is the standard plaquette action: 

Sg = ~P ]T ReC/ 12(s), (14) 

s 

with Uxz{s) = Ui(s)U 2 {s + l)U\{s + 2)U ] 2 {s). In the unquenched case the violation 
of LITh can certainly not be worse, since the effect of the fermion determinant will 
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Figure 2: Distribution of the values of Nr — Nl (full circles) and N r /2 d (empty squares) in 
a lattice of L = 8 at (3 = 1.0, compared to the distribution of the geometrical charge (solid 
line) . 

tend to suppress the measure of those configurations with non-trivial topology. For a 
similar study of the standard Wilson action in the unquenched case see, fl4|| . 

In the quenched Schwinger model, the distribution of the geometrical topological 
charge can be computed analytically. We recall the definition of the geometrical charge 
for this model: 

1 

2^ 



Q 9eo ^^J2 e ^( s ) = J2 lo s(Ui2(s)), \9 12 \<7T. (15) 



This quantity is an integer. The fraction of quenched configurations with a given 
geometrical charge Q is then given by 



Z 2ir J-tt I Q (u) ttv + a/2 

where Z = Iq((3) l2 and I u (x) are the usual modified Bessel functions and L is the 
number of lattice sites in each direction. Of course, we have checked that this dis- 
tribution is correctly reproduced by our Montecarlo routine. In Fig. 2, we compare 
this distribution to that of Nr — and N r /2 d calculated numerically for /3 = 1.0 
and 2 x 10 3 configurations (a ~ 0.3 in physical units, with respect to the mass gap in 
the continuum, i.e. m = l/t/W/3). The quantity N r /2 d , where iV 1 " is the total num- 
ber of real eigenvalues, has been proposed by some authors as a possible definition of 
topological charge |14}| . Clearly at this value of j3, it is very distorted near zero. 
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1 


8 


0.588 


0.186 


1.596 


0.281 


1.6 


10 


0.484 


0.192 


1.701 


0.261 


2.3 


12 


0.364 


0.191 


1.818 


0.216 


3.0 


14 


0.242 


0.133 


1.902 


0.144 


4.0 


16 


0.151 


0.057 


1.949 


0.074 


6.0 


20 


0.091 


0.022 


1.981 


0.017 



Table 1: Real Eigenvalue Distribution in lattice units for different (3. 



On the other hand, the distribution of Nr — Nl is also narrower than the geometrical 
one. This could be interpreted as due to the fact that the geometrical definition can 
measure the topological charge of smaller lumps of charge, being then more sensitive 
to dislocations. Let us suppose that the correct charge distribution was that obtained 
from the counting of real fermion modes. For each given configuration, the geometrical 
charge will differ from the fermionic one by the presence of one (to leading order, 
since a larger difference would be more rare) small lump of unit charge of either sign. 
This means that the configuration will be counted with equal probability as having 
a geometrical charge smaller/larger in one unit with respect to the fermionic charge. 
Since the charge distribution is always peaked at zero and decreases for increasing 
charge, then more configurations get shifted outwards than inwards, broadening it. 

However, this broadening could also be explained in terms of the mixing between 
the physical and doubler sectors of real eigenvalues. Let us now suppose that the good 
distribution is the geometrical one; then according to the LITh, the net chiral charge 
in S p is Q geo . If Op^o^ are large (which will happen at strong coupling) the modes 
corresponding to the first doubler sector can get inside the region defined as S p or 
vice versa. Taking into account that the net chirality in the first doubler region is 
opposite to that in the physical region, it is easy to convince oneself, by considering 
several examples, that the leading effect of this mixing is to narrow the fermionic charge 
distribution with respect to the geometrical one. If this is the mechanism by which the 
fermionic charge is not sensitive to very small lumps of topological charge, it is clearly 
not very justified to say that the fermionic charge has less lattice artefacts. 

As we move towards smaller coupling, the three charge distributions become much 
closer, although there is always a systematic difference between N r /2 d and the other 
two. As we have argued, in the Schwinger model, when (3 is increased, the chiral 
properties of action improve like 0(a/(3). In Table 1, we present, for different [3 values 
and fixed physical volume, the average and dispersion of the real eigenvalues (in lattice 
units) in the physical and first doubler regions. There m c = (A) p , A — 2r = (A)^, 
and a v fall as 1 / [3 for small enough coupling, confirming the expectation that they are 
0(a) effects. Notice that the values of a p ^ are misleading at large couplings, because 
the corresponding distributions are far from being Gaussian (see Fig. 1). 
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Figure 3: Probability to find LITh satisfied as a function of at fixed physical volume, 
for three values of the physical volume. 

In Fig. 3 we present the probability to find the LITh being satisfied at fixed physical 
volume, for three values of this volume. Each point corresponds to 200 configurations 
(separated by 100 Monte Carlo sweeps for the smaller (3 values and by 500 for the 
large (3 values). The probability grows for smaller as shown in Fig. 3 up to a 
point 1/Pc where it saturates to 1. The important point is that (3 C does not seem to 
depend on the physical volume. The quantities a p and (X) p also vary very little with 
the physical volume. An interpretation of Fig. 3 is that there is some critical splitting 
of the doubler sectors which ensures the validity of the LITh. 

If the violations of LITh were related to small objects such as dislocations, one would 
expect a dependence of (3 C on the volume, since the entropy and thus the Boltzman 
weight of these artefacts grows with the volume at fixed (5. In this model there are no 
dangerous artefacts at weak coupling (e.g. dislocations), because there are no scale- 
invariant instanton solutions and, as the continuum limit is approached, small lumps 
of topological charge are expected to be strongly suppressed. So if we want to draw 
any conclusion about the effect of dislocations in four dimensions on the LITh from the 
results in this two-dimensional model, the action should be improved while remaining 
at arbitrary strong coupling. This is possible with the action discussed in the next 
section. 
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4 LITh for the Improved Action 



From now on, the lattice spacing is denoted by b to distinguish it from the lattice 
spacing on the fermion lattice / and from the lattice spacing for the standard action 
a. The /-lattice is some integer subdivision of the 6-lattice. The b sites are denoted by 
s, while the sites on the finely-grained lattice are x. The path integral is defined as 

Z[rj\ = J VU^ e~ 5 ^] det(P) e"^" 1 ", (17) 

where rj are external fermion sources. The gauge action is the standard plaquette 
action of eq. (|i~4| ) and the measure is the standard one in a lattice of spacing b. The 
only difference is in the Wilson-Dirac operator, which is defined on the /-lattice as 
follows: 

P - lj2l,H + D-)- r -j2D+D-, (18) 

where the covariant and normal derivatives are given by D^^f(x) = u^(x)^>(x + ft) — 
^(x), D~^(x) = i S(x) — uj l (x — jl)'^f(x — p,). The u^x) link variables are interpolations 
of the real dynamical fields L^(s) ||. The reader is referred to for details on the 
method to construct the interpolation for non-Abelian theories. In this U(l) model 
the interpolation is particularly easy. Defining the sites on the /-lattice as x = sb + 
ix/1 + £2/2, where < t\,t 2 < 1, the interpolation that we use for compact U(l) in 
two dimensions is given by, 

/ f 

Ui(ti,t2) = J(h,t2) Pl 2 (s)] t2 b W(ii + ^,t 2 ) 



f 

U2{t\,t 2 ) = <J(t 1 ,t 2 ) uj(h,t 2 + 



(19) 



where the u> fields are defined as 



w(0,t a ) = U 2 (s)^, 

u(l,t 2 ) = C7i(s) (C/ 2 (s + i)) t2 , 

w(tl,*2) = (UiisY 1 ) 1 -' 2 [U 12 (s) h U 2 {s)U l {s + 2t\ t \ (20) 

It is trivial to show that this interpolation is gauge covariant, rotationally invariant up 
to a gauge transformation, and strictly loca]0. 

The continuum limit of this theory is defined as 

r X b^0 (21) 



2 Notice that this interpolation is not the same one proposed in q|. The reason is that ( fL9[ ) is not 
appropiate for chiral gauge theories, because the interpolated gauge held is not bounded. 
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for f/b fixed. The real cutoff is b. We can actually construct the improved fermion 
operator in terms of U^, taking any number for f /b, without ever referring to the 
lattice /. 

From the power counting arguments of J?J, the breaking of chiral symmetry is 
suppressed by at least a power of the ratio f/b to all orders in the gauge coupling. In 
the appendix, we compute the one-loop correction to the quark mass renormalization 
in QED in four dimensions, to leading order in f/b, and show this suppression by 
explicit calculation. It is not hard to understand this suppression along the lines of 
section 2. The continuum field can easily be chosen to also satisfy = exp(i / a M ). 
The operator ( |l8| ) can be expanded in / to give 

P=W + f D*D* + ... =P C + 0(f/b) (22) 

The higher-dimensional operators are truly of O(f), in contrast with the standard 
action, because the roughness of the continuum field is at most of 0(1/6). Similarly, 
close to the doubler momenta, 

= r + f D-D* + ... = -p+ C + 0(f/b). (23) 

For finite f/b, (P4 ) implies that the LITh is satisfied for arbitrarily large coupling, 
by the arguments of section 2. Furthermore, the real modes in S p have magnitudes 
of 0(f /b), without any tuning to the chiral limit and their quiralities are ±1 up to 
0(f /b) corrections. 

All these expectations are confirmed numerically. In Fig. 4 we show the dispersion 
of the real eigenvalues in S p , for a simple fine-graining factor of 1/2. The values have 
been shifted to an average value of zero. The improvement in the magnitude and 
dispersion of the real eigenvalues is roughly (f/b) 2 (in / units), as expected. This is in 
contrast with the results of ||, concerning the Clover improved action for this model. 

On the other hand, the improvement on the LITh is much more dramatic. In Fig. 5, 
we show the probability of finding the LITh, i.e. Q geo = Nr — Nl, as a function of 
(3 for a fr-lattice Lf, = 8 and a fine-graining factor of 1 (i.e. f = b = a) and 1/2 (i.e. 
/ = b/2). The number of uncorrelated configurations is 500 (separated by 500 sweeps) 
in each case. No violations of the LITh have been observed for the improved action 
for arbitrary rough fields. Also for fixed (3, the probability remains 1 at larger physical 
volumes. On the other hand, it is worth pointing out that, even for the improved 
action, the charge defined as N r /2 d does not coincide with the geometrical one. For 
example, the probability of their agreement is 0.93 at (3 = 0.01, improving slowly for 
larger (3. 

The fact that the LITh is satisfied for the improved action at arbitrarily strong 
coupling indicates that, as expected, similar results should hold in four dimensions. If 
this is so, there is no reason to believe that the fermion method to measure topological 
charge has less artefacts than the geometrical one. The problem of dislocations is 
related to the gauge action, which is very poor at measuring the action of small sized 
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-0.5 -0.4 -0.3 -0.2 -0.1 



Figure 4: Distribution of real eigenvalues (in / units) for f/b = 1/2 (full line), f/b — 1 
(dashed line) for a lattice size of = 8 and for (3 = 1.0. 



objects carrying topological charge, and the sensible way to get rid of this problem is to 
improve the gauge action [11, 12, 13]. On the other hand, measuring topology through 
the counting of the small real eigenvalues (taking into account their chiralities) might 
be computationally more efficient. A study of this issue will be presented elsewhere. 



5 Discussion and Conclusions 

It is well known that the quenched approximation has problems near the chiral limit. 
In the so-called exceptional configurations, it is very hard or impossible to invert the 
fermion matrix. Recently, a cure has been proposed for this 0]. Exceptional config- 
urations occur whenever there are real eigenvalues that get close to zero. As is clear 
from Fig. 4, this is bound to happen if the subtracted fermion mass is smaller than 
dp. The proposal of || is that whenever an exceptional configuration is found, the 
real eigenvalues must be shifted by the appropriate amount to be exactly zero at the 
chiral point, which implies that the fermion matrix is invertible arbitrarily close to this 
point. As already stated in ||, this is a non-local procedure, which might introduce 
unphysical effects. It is unclear how one could be sure about the safety of this method, 
even if it seems to give good results for a given observable. 

In the light of our results, it is important to stress that exceptional configurations 
are simply topologically non-trivial configurations. Clearly the quenched approxima- 
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■ f/b = 1 /2 



D.5 1 1.5 2 2.5 3 

S 

Figure 5: Probability of finding Q 9eo = Nr — Nl in a lattice of Lb = 8 as a function of (3 for 
f /b = 1 (circles) and 1/2 (squares). 

tion should not be well-behaved in non-trivial configurations, since a big suppression 
in the Boltzman weight of such configurations comes from the fermion determinant, 
which is being neglected. A much more justified procedure would be to ignore topo- 
logically non-trivial configurations in this approximation. Alternatively one should not 
approach the chiral limit closer than a p (0(f /b) 2 for the improved action or supposedly 
0(a) for the standard action), i.e. work with a quark mass > a p . This is a sensible 
thing to do even for the standard action if, as argued in |J the dispersion of the 
real eigenvalues is a lattice artifact that vanishes in the continuum limit. In other 
words, the physical fermion mass cannot be determined with a better accuracy than 
0(a) |flO| ], so working with a physical mass of order a p = 0(a) is not surprisingly the 
safe choice. 

The lesson is that the physical mass should then not be reduced without improve- 
ment. At least in principle f/b can be taken to be as small as needed []. Our improved 
action then has the same effect as the "modified" quenched approximation, i.e. one can 
get closer to the chiral point without encountering exceptional configurations, except 
that it does it in a local way. For another possible solution to this problem, see |16[| . 

For full dynamical fermions, the problem of exceptional configurations is not so 
important, because they have a very small Boltzman weight. Of course this means that 

3 Of course lowering f/b is hard computationally, but it should be possible to parametrize the action 
for arbitrarily small f/b in terms of a b lattice fermionic action with higher dimensional operators with 
coveniently tuned coefficients. 
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they are less probable, but can and will appear. On the other hand, it is clear that 
configurations with real eigenvalues are a serious problem for unquenched simulations 
with an odd number of flavours. Figure 4 shows that if the renormalized mass is smaller 
than 0(a p ), we will often find configurations with a real negative eigenvalue (notice 
that complex eigenvalues come in complex conjugate pairs, so their contribution to the 
determinant is always positive), which makes the determinant negative; the Euclidean 
measure is thus not positive-definite and Montecarlo methods will fail. Again the 
obvious solution is to choose a quark mass of 0(a p ) or treat the negative eigenvalues 
as in [16]. 

To conclude, we have presented evidence for the exact validity of the index theorem 
on the lattice if the zero modes are identified with the exactly real eigenvalues of an 
improved Wilson-Dirac operator in the physical region (i.e. A < r/a) and if the gauge 
topology is measured with the geometrical definition of ref. [Qj. We have argued that 
the violations of this "theorem" that are observed for rough configurations are related 
to the large dispersion of the real eigenvalues: the different doubler regions cannot 
be cleanly separated in the real axis and the counting of the physical chiral modes 
gets contaminated by the doublers. The effect of improving the fermionic action is to 
ensure a finite gap between the doubler sectors, which then ensures the exact validity 
of the "index theorem" for arbitrarily large coupling. These results strongly suggest 
that fermionic methods of measuring an integer-valued topological charge using the 
Wilson operator are equivalent to the geometrical definition. 

A few days ago two new preprints appeared which study the index theorem for 
two other types of improvement: perfect actions [17] and Clover improvement in the 
second Actually our improved fermionic action in the limit f /b — > 0, combined 
with the improved gauge action of [13], shares many properties with perfect actions. 
Not surprisingly the index theorem becomes exact in both cases at arbitrarily large 
gauge coupling. 
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6 Appendix A 

In this appendix, we sketch the one-loop calculation of the mass renormalization in 
QED in four dimensions with the two-cutoff action. In perturbation theory, the in- 
terpolation of H simplifies greatly, and an analytic formula can be obtained, which 
corresponds to a linear interpolation^. The Fourier transform of the interpolated in 

4 Notice that in lattice perturbation theory there are no problems with the windings discussed in | 
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I 1 # 



Figure 6: Diagrams contributing to the fermion self-energy at one loop. 



field in four dimensions is, 

1 , N e^ n -1 1 e iqan 2(cos z a - 1) . 

<*<«> = ^M-Strr 5 ^ (e i- iP < 24 > 

where q u = q~u + ^ h± , \%\ < % and the integers m u = -n/2, ...,n/2. z M = <? M n = c^n 
Notice that at low momentum, a M (g) — > ^(<?), as expected from the locality of the 
interpolation. 

In order to derive perturbation theory, we need to go to momentum space. Except 
for the gauge-field propagator, all the lattice Feynman rules are the same as for standard 
Wilson fermions on the / lattice. We will use / = 1 units to simplify the formulae. 
Then we should extract the leading dependence on n = b/f, which is the factor of 
fine- graining. 

From (p4[), the propagator of the interpolated gauge field can be obtained in terms 
of the gauge propagator on the b lattice in the Feynman gauge: 

MgM-g)> = 5, " 2 JJ-^gffi J] i(i-^) a . (25) 

2 2^p sm(z p /2)^ nP 1 — cos q a V" n z 1 — cos q a 

Now, we want to compute the linear divergence in the fermion self-energy. The two 
contributing diagrams are those shown in Fig. 5. In order to get the linear divergence, 
the external fermion momenta can be set to zero. The tadpole contribution is given by 

2 1 f a z n s—^ 1 ^ 1 — cos z a 



" 1 1 2 n 2 



S£l(0) ^ J BZ (2tt)4 2 E P sin(z p /2)2 E n« 1 " cos ^ 

n/2 

n e ±c-^^r. (26) 

, /n n z 1 - cos q a 

The integer sums can be analytically computed, using 

n/2 

E 



m M =— n/2 



sin(g M /2 + irm^/n) 2 1-cosz^' 



™/2 -, A r, , 

y 1 _ xl 4 4 2 + cosz M 

v _ n/2 sin (^/2 + vrm A1 /n) 4 3 (1 - cos z^) 2 ' 
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Finally, we get 

S (0) " ~ 9 Wrf Lz \ 2E p sin(z p /2)2 + O(^). (28) 

The remaining integral is a finite number, which can be computed numerically. The 
leading dependence on n is thus the expected one. 

The contribution from the second diagram is a little harder to obtain, because the 
fermion propagator enters the integer sums; however, the leading dependence on n can 
also be obtained analytically. The diagram gives 

y b (()) = -n 2 — ( d ^ z n2 V 1 V 1 ~ C0S ^ 

1 ] ~ 9 n*J BZ (2tt)4 2 E P sin( z p /2)2 ^ n * m ^ n/2 1 " cos % 

n 1 1 - coszq 2 cos q p M(q) - sin 2 g M 
, { ^ ,n 2 4-cosg a j ' s(q) + M(qy ' 1 ^ 

a^tfi m a =—n/2 ^ w/ 

with M(q) = J2 P 1 ~~ cosq p and s(q) = J2 P i s ^ n Qp) 2 ■ The integer sum over m p can be 
derived from eqs. (p7|) after some manipulations: 

v cosg^M^-sin^ _ 2 1 M'(g) 

^ (1 - cos q,,) (s(q) + M(g) 2 ) 1 - cos s'(q) + M'(g) 2 lJ lJ 

where s'(q) = J2p^p( s ^ n Qp) 2 an d M'(q) = J2 P ^p 1 ~~ cosq p . The next sum required is 
V M'(q) 1 n 4 2 + cosz. M"(g) 

^ s'(g) + M'(g) 2 (1 - cos g^) 2 3(1- cos z v f s"(g) + M"(g) 2 + W ' 1 J 

with s"(q) = J2p^^,u( s ^ n Qp) 2 an d similarly for M"{q). The last two sums are also of 
the form ([31]) and the final result is 

*®--*i±?L£? ^^ n (2+a * !j+01 ^' (32) 

Again the remaining integral is a finite number, so the n dependence is as expected. 
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